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REMARKS 

Claims 10-14, 16-20, 25 and 27-50 were pending in the present application. Claims 11, 
25, 30, 33, 38 and 45 have been cancelled without prejudice. New claims 51-54 have been added. 
The new claims are supported in the claims as originally filed, and at least in the specification at 
page 8, lines 24-27. Thus, claims 10, 12-14, 16-20, 27-29, 31-32, 34-37, 39-44 and 46-54 are 
pending. 

To further prosecution without acquiescing to the Examiner's arguments, Applicants 
have amended independent claim 10. As amended, claim 10 relates to a process for preparing 
randomly-ordered crystal agglomerates comprising an alkali metal clavulanate salt, comprising 
contacting a solution or suspension of alkali metal clavulanate salt in a solvent or mixture of 
solvents with one or more anti-solvents under stirring. This amendment is supported in the 
specification and does not include new matter. (See e.g., specification at page 5, lines 21-28; and at 
page 7, line 30 through page 8, line 2). Applicants request reconsideration of the Examiner's 
rejections in the Final Office Action mailed March 10, 2004, maintained in the Advisory Action 
mailed April 6, 2004, in view of the amended claims. 

Rejections under 35 U.S.C. $ 102 

First, the Examiner rejected claims 37-38, 40-42, and 44-50 under 35 U.S.C. § 102(b), as 
allegedly being anticipated by U.S. patents 4,454,069; 6,417,352; and 5,288,861. Applicants 
address the Examiner's rejection in view of the amended claims. As amended, claim 37 relates to 
an agglomerate of randomly-ordered crystals of an alkali metal clavulanate having a compressibility 
between about 10 % and 40 %, with the proviso that the rosette-like crystalline form of potassium 
clavulanate is excluded. Claims 37, 39-44, and 46-50 are not anticipated because U.S. patents 
4,454,069; 6,417,352; and 5,288,861 only describe well-defined needles or plates, or rosette forms. 

U.S. 5,288,861 ('861) describes potassium clavulanate in rosette form, where a plurality 
of needle crystals radiate out from a common nucleation point. The '861 patent also teaches that 
crystalline potassium clavulanate generally exists in the form of rod-like or needle-like crystals. 



sd-l97990vl 



Application No.: 09/937,834 



7 



Docket No.: 246152015300 



The crystals are sometimes agglomerated into plate-like crystals, or sometimes randomly 
aggregated into' loosely formed bundles. These forms of potassium clavulanate can give rise to 
processing difficulties because the material is of low bulk density, and does not always flow readily. 
('861, at col. 1, lines 40-48). Contrary to randomly-ordered crystal aggregates, needle crystals by 
definition, are well-defined structures having a general form. (See e.g., Langer, "Existence of 
needle crystals in local models of solidification, Phys. Rev, A 33: 435-441, attached as Exh. 1). 
Furthermore, as rosette-like crystalline forms of potassium clavulanates have been excluded in 
claim 37, this patent does not anticipate the claimed agglomerates. 

U.S. patent 4,454,069 ( c 069) teaches a process for preparing potassium clavulanates as 
well-defined needles or waisted plates (i.e., butterfly-shaped). ('069, at col. 5, lines 32-39). The 
clavulanates in the '069 patent were prepared by adding potassium 2-ethyl hexanoate to a solution 
or suspension of an amine-clavulanate salt. 

U.S. patent 6,417,352 ('352) teaches a process for preparing potassium clavulanates that 
is similar to the '069 patent, where potassium-2-ethyl hexanoate is added to clavulanic acid, instead 
of an amine-clavulanate salt. As indicated at column 4, lines 40-56, "conversion from sodium 
clavulanate to potassium clavulanate is carried out by extraction of clavulanic acid to an adequate 
solvent and crystallization of potassium clavulanate after dilution of the acid . . . and addition of 
potassium 2-ethylhexanoate or potassium acetate solution." Although the '352 patent does not 
specify the form of potassium clavulanates, conditions which would be expected to obtain needles 
forms were used Because the * 861, fc 069 or *3S2 patent do not teach aggregates of randomly- 
ordered crystals of an alkali metal clavulanate salt, claims 37, 39-44, and 46-50 are not anticipated. 
Applicants therefore, respectfully request that this rejection be withdrawn. 

Second, the Examiner rejected claims 37-50 under 35 U.S.C. § 102(b), as allegedly 
being anticipated by WO 97/33564. Claims 37, 39-44, and 46-50 are not anticipated because WO 
97/33564 only describe agglomerates of penicillin V potassium, phenoxymethylpenicillin 
potassium, amoxicillin trihydrate, and cephalexin monohydrate. While these antibiotic 
agglomerates may be mixed with a second pharmaceutical^ active agent such as potassium 
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clavulanate, the potassium clavulanate itself is not an agglomerate of randomly-ordered crystals. 
Thus, WO 97/33564 does not anticipate claims 37, 39-44 and 46-50, and Applicants respectfully 
request that this rejection be withdrawn. 

Third, the Examiner rejected claims 10-14, 16-19, 27, 27-32, 37-38 and 42-45 under 35 
U.S.C. § 102(b), as allegedly being anticipated by WO 98/21212. Claims 37, 39-44, and 46-50 are 
not anticipated because WO 98/21212 only describes rosette or needle formed crystals, or a cluster 
of needle formed crystals (WO 98/21212, at page 4, lines 17-18). As previously indicated, needles 
have well-defined structures, unlike the randomly-ordered crystal clusters of the claimed invention. 
Thus, WO 98/21212 does not anticipate claims 37, 39-44 and 46-50, and Applicants respectfully 
request that this rejection be withdrawn. 

Furthermore, the needles obtained using the process taught in WO 98/21212 have a 
compressibility of about 50 %, as shown in Table 1 (Specification at Example 8). In the Advisory 
Action, the Examiner indicated that w [t]here is no way of knowing where this needle material came 
from." (Advisory Action, page 3). Contrary to the Examiner's assertions, the needles were 
prepared using conditions that are similar to those described in Example 6 of WO 98/21212. In 
particular, needles of potassium clavulanate were prepared by suspending diclavulanate salt of 
bis(2-dimethylaminoethyl) ether in acetone and water. Under stirring, a solution of potassium 2- 
ethylhexanoate in acetone was added. (See, specification at page 1 7, lines 4-1 0). 

WO 98/21212 also fails to teach a process for preparing crystallized agglomerates 
comprising an alkali metal clavulanate salt, comprising contacting an alkali metal clavulanate salt in 
a solvent or mixture of solvents with one or more anti-solvents under stirring. In contrast, WO 
98/21212 teaches a process where a potassium source in a solvent, is added to a amine salt in a 
solvent, preferably in a water miscible ketone or alcohol. In a particular example, potassium 2- 
ethylhexanoate in acetone was added to bis(2-(dimethylamino)ethyl) ether diclavulanate in a 
mixture of acetone and water. In the Advisory Action, the Examiner indicated that "[t]he exact 
instant the K salt is formed in that solution, it will meet the claim language requirement of 
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potassium clavulanate, and at that point, the solvent limitation is met as well." (Advisory Action, 
pages 2-3). 

Contrary to the Examiner's assertions, the formation of potassium salt in solution alone 
does not meet the claim language, because the potassium clavulanate solution is not contacted with 
another portion of anti-solvent. On the other hand, agglomerates of randomly-ordered clusters of 
alkali metal clavulanate salts may be prepared, for example, by adding a solution of potassium 
clavulanate in a mixture of water/acetone to acetone. (See e.g., specification at Example 1). 
Furthermore, as previously indicated, the potassium clavulanates obtained using the process taught 
in WO 98/21212 are rosette or needle formed crystals, or a cluster of needle formed crystals. 
Because WO 98/21212 fails to teach agglomerates of randomly-ordered clusters of an alkali metal 
clavulanate salt, or a process of preparing such agglomerates, claims 10, 12-14, 16-20, 27-29, 31-32, 
34-37, 39-44 and 46-50 are not anticipated. Applicants therefore, respectfully request that this 
rejection be withdrawn. 

Further, new claims 5 1 -54 are not anticipated by the prior art references cited by the 
Examiner. Specifically, claims 5 1-52 are dependent on claim 37, and contain all the limitations in 
claim 37. Claim 53-54 are dependent on claim 10, and contain all the limitations in claim 10, As 
previously discussed, claims 10 or 37 are not anticipated. Accordingly, claims 51-54 are not 
anticipated. 

Rejection under 35 U.S.C. S 1 12, second paragraph 

The Examiner rejected claims 31 and 45 under 35 U.S.C. § 1 12, as allegedly being 
indefinite, for reciting "compressibility" (claim 3 1) and "flowability" (claim 45). Claim 45 has 
been canceled without prejudice, rendering this rejection moot. 

Regarding "compressibility," Applicants have amended the claims to recite "Carr index 
compressibility." The recitation of the term "Carr index" merely includes art-recognized definitions 
known at the time of filing the application, and is not considered new matter. See MPEP § 2163.07. 
As commonly known in the art, Carr's index is the percentage calculated as 100 times the ratio of 
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the difference between tapped bulk density and loose bulk density to the tapped bulk density. (See 
e.g., Gohel et al. 9 "Development and evaluation of a multifunctional directly compressible diluent 
consisting of brittle and ductile materials," Pharm. Tech., December 4, 2003, attached as Exh. 2). 
The Can* index was also used for calculating compressibility in the Examples. As amended, the 
claims are definite, and Applicants respectfully request that this rejection be withdrawn. 

Rejections under 35 U.S.C. §112. fast paragraph 

The Examiner rejected claims 10-14, 16-20, 25, 27-32, 35-45 and 47-50 under 35 U.S.C. 
§ 1 1 2, first paragraph, as allegedly failing to comply with the enablement requirement. As 
amended, the claims relate to agglomerates comprising an alkali metal clavulanate salt, and 
processes of preparing such agglomerates. In the Advisory Action, the Examiner indicated that 
such an amendment, if entered, would result in enabled claims. Thus, Applicants respectfully 
request that this rejection be withdrawn. 

Furthermore, the Examiner rejected claims 10-14, 16-20, 25 and 27-50 under 35 U.S.C. 
§ 1 12, first paragraph, as allegedly failing to comply with the written description requirement The 
Examiner alleged that the removal of a previous limitation of "high water affinity" results in 
compounds beyond what the specification teaches. (Final Office Action, page 6). The Examiner 
also alleged that the removal of the previous limitation of "high water affinity" has broadened the 
claims. Applicants must respectfully disagree. While the previous limitation of "high water 
affinity" has been removed, the amended claims relate specifically to agglomerates of randomly- 
ordered clusters of an alkali metal clavulanate salt, and processes of preparing such compounds. 
The amended claims reasonably convey to one skilled in the art that the inventors) had possession 
of the claimed invention at the time the application was filed. Thus, Applicants respectfully request 
that this rejection be withdrawn. 

In view of the above, each of the presently pending claims in this application is believed 
to be in immediate condition for allowance. Accordingly, the Examiner is respectfully requested to 
withdraw the outstanding rejection of the claims and to pass this application to issue. If it is 
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determined that a telephone conference would expedite the prosecution of this application, the 



Examiner is invited to telephone the undersigned at the number given below. 

In the event the U.S. Patent and Trademark office determines that an extension and/or 
other relief is required, applicant petitions for any required relief including extensions of time and 
authorizes the Commissioner to charge the cost of such petitions and/or other fees due in connection 
with the filing of this document to Deposit Account No. 03-1952 referencing docket no. 
24615201 5000. However, the Commissioner is not authorized to charge the cost of the issue fee to 
the Deposit Account. 

Dated: May 6, 2004 Respectfully submitted, 
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Existence of needle crystals in local models of solidification 
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L INTRODUCTION 

There is growing evidence that pattern selection in local 
models of dendritic solidification is closely connected 
with the existence of so-called "needle-crystal" solutions 
of the steady-state equations of motion. In both the 
geometrical 1 ' 2 and rjoimdary-iayer* 4 models that have 
been introduced recently, as well as in more realistic non- 
local models of solidification with diffusion control, 5 
there exists a continuous family of "Ivantsov'*- 7 needle 
crystals when surface tension is omitted The addition of 
surface tension, however, is a singular perturbation of 
these systems which, at least in the local versions— that is, 
the geometrical and boundary-layer models—destroys the 
family of solutions and opens the possibility that sharp 
selection of growth rates, tip radii, etc., occurs via a lolva- 
bility condition. It seems possible that a similar 
phenomenon occurs in the fully nonlocal problem. 

The purpose of the present paper is to examine in some 
detail the breakdown of the Ivantsov solutions in both of 
the local models. The way in which this breakdown 
occurs provides some interesting clues about the roles 
played by surface tension and crystalline anisotropy in the 
dendrite theory. Moreover, this singular perturbation 
problem is of mathematical interest in its own right. As 
we shall see, the method of solution proposed here seems 
reasonable but is not rigorous or even systematic. It can, 
however, be tested by direct numerical computations; and 
the results of such tests make it plausible that the analytic 
approximation captures, the essential features of the 
relevant phenomena. Indeed, as this paper is being writ- 
ten, it appears that there may be important new progress 
both in understanding the mathematical nature of the 
problem as posed here' and in showing that the basic 
features of this model problem also appear in the more 
realistic systems. 9 -" This paper has therefore been or- 
ganized to serve as a starting point for study of these later 
developments* 

We shall begin in Sec II by looking at the simplest non- 
trivial version of the geometrical model in order to 
describe the mathematical strategy with a minimum of 



unnecessary complication. That strategy will be outlined 
in Sec. HI. In Sec. IV we shall apply this method to a 
minimal version of the boundary-layer model, and shall 
see there how crystalline anisotropy may play an essential 
role in the selection of needle-crystal solutions. Some 
mathematical and computational details pertaining to 
both tnddels are relegated to an appendix. 

H. THE QUASILINEAR GEOMETRICAL MODEL 

In all of the following, we shall restrict our attention to 
two-dimensional local models in which a moving one- 
dimensional interface is described by specifying its curva- 
ture K^dS/ds as a function of arc length 5. The 
relevant geometry is illustrated in Fig. L The basic as- 
sumption of the geometrical models 1 ' 3 is that v nt the nor- 
mal velocity of the interface, is a function only of K and 
its even derivatives with respect to s. The steady-state 
condition is simply 

o M iK t d 2 K/ds 2 f . . .)=t?cos0 , (2.1) 

so that the interface is moving at constant velocity u. 
without change in its shape, in a fixed direction. A needle 
crystal, by definition, is a solution of (2.1) that has the 
general form shown in Fig. 1 in which 




FIG. 1. Geometry of the needle crystal. 
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The analog of the Ivantsov limit in the geometrical 
model occurs when we set v 9 to be a function of K alone. 
If (2.1) is ptcccwise invertible for A" as a function of 0, 
then Ki9)=md9/ds is a first-order ordinary differential 
equation which generally is solvable for 0 as a function of 
s. We may further expect that these solutions will satisfy 
the needle-crystal conditions (2.2) for some continuous 
range of values of the parameter u. For example, the sim- 
plest possible choice is v n **K t for which it turns out that 

0=cos- , (sechw), (2.3) 

a relation which satisfies (2.2) for all u. The family of 
Ivantsov solutions is conventionally characterized by the 
relation between v and the tip radius R = AT" 1 (0=0), 
which in this case is trivially Rv ~ 1. 

The above solutions are all strongly unstable. In fact, 
the model as it stands is not even dynamically well de-' 
fined because the amplification rate of deformations is un- 
bounded at short wavelengths. In order to produce a 
model whose time dependence is meaningful, it is neces- 
sary to add a term which mimics the effect of surface ten- 
sion and stabilizes the system at short wavelengths. The 
easiest way to do this is to write 



ds 1 



(2.4) 



where y plays the role of a capillary length. Equation 
(2,4) defines what might be called a quasilinear geometri- 
cal model; the right-hand side is linear in K. This model 
does not actually produce dendrites, but its dynamical 
behavior is well defined and quite interesting in some 
respects. 

Our aim now is to see what effect the second derivative 
in (2.4} has on the needle-crystal solutions (2.3). To do 
this, it is convenient first to rewrite the steady-state equa- 
tion (2.1) in the form 



K=COS0— y 



d 1 * 



(2.5) 



(2.6) 



where v=(yv) l f K «=wr. and §=vs. Equivalently, 
2 d& 

The latter form, in which £ has been replaced by 0 as the 
independent variable by using *=d$/dg, is very useful 
but can produce spurious 4ifficuities unless k is every- 
where non-negative. 

Equations (2.5) and (2.6) can be used to illustrate some 
points that have been made in previous papers. 2 ' 3 If we 
iterate the right-hand side of (2.6), we generate a series ex- 
pansion for *r in powers of v and cos0. Each term in this 
series is consistent with the needle-crystal conditions (2.2), 
and the series appears to produce an accurate estimate for 
H especially far down the needle where cos0 is small. 
However, we know that this series can be at best asymp- 
totic. To see this, write (2.5) as a set of three coupled 
equations for a f -dependent trajectory in the space of vari- 
ables 0, K**d9/di, and As has been noted 
previously, 1 only one trajectory emerges from the fixed 



point at 0=* -ir/2, *<=a=0; similarly, only one trajectory 
enters the reflected fixed point at 0=+ir/2. For a needle 
crystal to exist, these two pieces must belong to a single 
trajectory which joins the fixed points and, by symmetry, 
passes through some point on the k axis with 0=A,=o' 
There is no special reason for this to happen* In general 
the trajectories entering or leaving the fixed points, even if 
well approximated by some finite number of terms of the 
above series in powers of v, need not reach 0»O— the tip 
of the needle— with A.«=</«/rf£«o. Our goal in what fol- 
lows is to compute k at 0=0 explicitly and to discover 
' under what conditions it might vanish. 

As we shall see immediately, the mathematical problem 
posed above is highly nontrivial To make some progress 
analytically, we shall look only at small values of v t that 
is, we shall look in the neighborhood of the known solu- 
tion of (2.6), *oW=cos0 at v«0. We then *h&ll consider 
solutions of (2.6) for nonzero v that satisfy the needle- 
crystal conditions (2.1) for 0— +tr/2, and shall examine 
their behavior near 6 =0. This strategy is based on the as- 
sumption that, if v is sufficiently small, the difference be- 
tween * and *q also will be small in the interval 
and may be computable by a linear approxi- 
mation. For arbitrary v, the continuation of this solution 
to negative 0 may not remain small or even well defined, 
but the latter behavior need not invalidate an approxima- 
tion in the region of interest. The trouble with this pro- 
cedure turns out to be that X(0=O) vanishes more rapidly, 
than any finite power of v. Indeed, a systematic expan- 
sion in powers of v must necessarily recover just the 
asymptotic series described above for which K(Q) vanishes 
identically term by term. As a result, any useful scheme 
of approximation must go beyond a simple series in 
powers of the small parameter v. 
. To examine this situation in greater detail, define 



#C— COS0 S VK\ , 



and write (2.6) in the form 



C.7) 



d 2 «i ^ d* x 
v^-2vtan0— +— n +v<l-3< 



dd 1 



1 



s 2 0)]k, 



1 

COS0 



U-2eort)+.^l*,), (2.8) 



where ^\k x ) is a nonlinear term of orders v 1 *? and v*k\ 
containing at most two differentiations of these quantities 
with respect to 6. The structure of (2.8) suggests that 
each differentiation with respect to 9 produces a factor 
1/Vv. We shall see this explicitly in what follows. Thus, 
a systematic approximation for (2*8) might involve drop* 
ping and the term proportional to v in square brackets 
on the left-hand side, or perhaps treating these terms per- 
turbatively. The term containing vd* s /d9, which is for- 
mally of order ✓v, turns out to be essential. In Sec III 
we shall examine the consequences of the above approxi- 
mation. It will be useful to describe the calculation in 
genera) terms so that the technique may be applied to the 
boundary-layer model in Sec. IV without repetitive ex- 
planations. 
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The general structure of (2.8), with or without the terms 
to be omitted, is 



de 



r{$) 



(3.1? 



The range of 9 is [0,tr/2]. Boundary conditions are to be 
imposed only at 0=ir/2. These are the needle-crystal 
conditions (2.2) which require that * and k 
= dK/d£=*Kd*/d$ vanish as £>+oo, 0—ir/2. If we are 
working only to fust order in *„ then this condition re- 
quires that *! and K 0 d* t /d9 vanish at ir/2, where * 0 
denotes the Ivantsov solution; that is, tf o =cos0 for the 
geometrical model Another way of stating this is that we 
want *| to be consistent near it V2 with the asymptotic ex- 
pansion of k in powers of v and ccsd. To lowest order in 
v, this means 

Ki**r 0 ($)/q 0 {9)> $-+ira 02) 

where r 0 and g 0 are the values of r and q at v=0. Final- 
ly, the problem that we pose for ourselves is to compute 
k=dK/dg 9 or simply dx/d9 9 at 0=0. 

For sufficiently small v t (3.1) can be solved by a 
WentzeJ-KramefS-BriUouin (WKB) approximation. The 
homogeneous solutions are 



9o 



tons are 



(3.3) 



(3.4) 



where 

and we have neglected terms of order Vv or smaller in the 
exponent and in the pref actor. These solutions should be 
accurate as long as 

1/2 

(3.5) 



1 dg 0 






g 0 d9 


« 


V 



farall9Jn[G>/2]. 

For the quasilinear geometrical model, $ o «=(cos0)~ 2 > 
and (3.5) is satisfied for v«l. Moreover, ijko) in (3.4) is 
equal to the unperturbed dimensioniess arc length £ ob- 
tained from * o =</0/d£=cos0. Using p 0 =— 2tan0, we 
find 



1 



i 



1-fsinfl 



cos0 



II 



'(coshf^exp j±^j . (3.6) 



For small v, these are rapidly oscillating functions whose 
amplitudes diverge as 0-nr/2, Equation (3.6) is 

an analytic statement of the fact that only one trajectory 
in $,K,k space enters the fixed point at 0=ir/2 f K=kr=0, 
and all other trajectories spiral out and away from that 
point. It also justifies our expectation that differentia- 
tions in (2.8) produce factors of 1/i/v, 
The next step is to use the homogeneous solutions (3.3) 



to compute a particular solution of <3.1).by means of the 
formula ' 



•-J^k^)], (3J) 



(3.8) 



where the Wronskian IV is 

^ d* ~df ' 

The lower limit of integration remains to be chosen. Be- 
cause both of the homogeneous solutions *$? ,n (0) which 
appear in (3.7) are inconsistent with the needle-crystal 
conditions, the only possible choice for this lower limit is 
ir/2. Using (3.3), we find 



(3.9) 



and 



Xexp(^4/^p o£ f^] 

Xsinf(l/v^)[#0)-#*)]} . Q.iO) 

To check that (3.10) does satisfy the asymptotic condition 
(3.2), integrate (3.10) twice by pans, integrating the sine as 
if to obtain a series in growing powers of v. The result is 



qoK9) 



(3,11) 



where M is a function of 6 which is of no special interest 
except for the fact that it exists and is well behaved near 

Our goal is to calculate dx/dO at 9—0, which we 
denote by the symbol *'(0). Note that the symmetry of 
the system requires that q and r be even functions of 9 
and that p be odd. Thus 



r o (0) 



x«p(t / 0 %<r>] 



Xcos 



(3.12) 



Equation (3.12) Is the principal result whose validity and 
implications are to be explored in the remainder of this 
paper. 

Remember that the condition for existence of a needle 
crystal is that k'{0) vanish. It is interesting, and possibly 
useful for future analysis, to note that (3.12) with k'(0)«0 
can be interpreted as a solvability condition for the (for- 
mally) linear inhomogeneous equation (3.1). To see this, 
consider (3.1) in the whole interval (— ir/2, +ir/2J and 
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is the WKB approriraation for the solution of ^ir«0 

5^ the linear operator on the left-hand side 

of (3.1), and 



1 



may 



V~y J 



(3.14) 

U tile solution of the adjoint equation ^=0. Thus, the 
solvability condition derived above has the form 

(3,15) 

wkfcfa is the usual statement that the inhomogeneous term 
r{9) can have no projection onto the null space of Jf if 
(3.1) is to have a solution. One obscure aspect of this in- 
terpretation is the definition of the function space: the 
divergent functions would not seem to belong to 

the space of acceptable needle-crystal solutions. One can 
make this difficulty seem less severe by working with the 
bounded functions that is, by transforming away 

the first derivative in <3.l>. But this mathematical point, 
among many such points raised in this paper, requires 
further scrutiny. 

The general character of (3.12) can be seen by evaluat- 
ing it for the quasilinear geometrical model where 
p o = -2tan0, r o *=(2co8 2 0- U/cos0, 
^Zh^ c «0=s«h£. Then, transforming to £ as the 
variable of integration, we find 



*'(0)c*-! *£<sech 3 ^)(2sech 2 £-l 



2Vv 



(3.16) 



where 



16%/^ 
15 



■^1.8906 



(3.17) 

The derivation of the second form of (3.16), valid in the 
lima v « 1, is summarized in the Appendix. The impor- 
tant point to notice is that (3.16) has an essential singular- 
ity at v~0, consistent with the failure of expansions in 
powers of v. 

A comparison between (3.16) and a direct numerical in- 
tegration of the fully nonlinear equation (2.5) confirms the 
basic form of (3.16), both the exponential function and the 
power of v in the prefactor. The numerically determined 
value of A, however, appears to be about 3,7±<U, larger 
than that given by (3.17) by a factor of about 2. this 
comparison was made for values of v in the interval 

? *Z? 0 0it acroM whicb *'< 0) charges by eight dc 
fades. The validity of the asymptotic evaluation of the 
integral in (3.16) was also checked in this interval by 
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direct numerical integration, thus confirming that these 
value* of v are small enough to be in the aayoX^ 

the T i^ r 5° f « related to the fact that 

the power of v in the prefactor in (3.16), i.e., v-»*u 
denned by the highest power of seel* in the in 
tegrand, which in turn can be traced back to (he highest 
power of cos* in , oW « O A0). (This will be seenTpU 

because there we are computing the high-frequency— 

about a first nonlinear correction to r. That is, consider 

»f, Wl M l tncIud f« result in r on the right-hand side 
? i fi' , ^ e J esuIhng orations to r will be formally of 
rdaUvely hjgh order in v but wiB contain varicT^e* 
«h. fX. ,atte S wi » P«««uce inverse powers ofv Z 

3 16)^.1^^" f ° f * m The "potential part of 
<3.16)wiH be unchanged by this procedure, but apparently 
theprefactor a not being computed systematically 

This one of analysis will not be pursued further in this 
Paper. It appears that there may be a much more elegant 
and mathematically controlled way of arriving at formu- 
la i such as (3.16). 8 The interested reader-along withTe 
of pr^s * * ,en 0CW devel °P raeats in this "las, 

IV. THE MINIMAL BOUNDARY-LAYER MODEL 

The boundary-layer model 3 - 4 has been introduced as an 
attempt to include in a local description some physical 
features of the real solidification problem that are missing 
in the purely geometrical approach. In particular, the 
dynamics of a thermal boundary-layer field defined along 
the solidification front mimic some of the nonlocality and 
history dependence associated with a more realistic 
thermal diffusion field. The model has had some 
successes. Its steady-state solutions, including the para- 
bolw needle crystal in the Ivantsov limit, are nearly identi- 
cal to those of the full model. The boundary-layer model 
does produce time-dependent dendritic patterns whose 
growth rates and tip radii are numerically consistent with 
solvability conditions for needle crystals with nonvanish- 
mg surface tension. At present, however, it is still not 
known whether the model produces a physically realistic 
picture of dendritic sidebrancbJng. 

In order to make the following analysis reasonably 
tractable, we shall consider only a minimal version of the 
boundary-layer model in which the dimenstonless under- 
cooling A becomes vanisbingry small. This is not the 
physically most realistic limit of the model, but it retains 
most of the features which seem essential for present pur- 
poses. (The opposite limit, A— »1, will be discussed else- 
where;) Th e steady-state equation to be considered is the 
same as that of Ref. (3), Eq. (5.17), supplemented by an 
anisotropic kinetic attachment coefficient. Specifically, 



«=cos*e— vxcos'd 



dB J cos 



d 



cos£ d$ 



(*+*«cos0) 



(4.1) 



33 



EXISTENCE OF NEEDLE CRYSTALS IN LOCAL MODELS OF . 



where the k in Ref 3 has been replaced by v* in order to 
emphasize Che analogy to (2.6). Here, v=d 0 u/D&*, where 
f** 3 » capillary length proportional to surface tension, D 
is the diffusion constant in the liquid, and v is the growth 
velocity as before. Note that, as in the geometrical model, 
v vanishes tn the limit of vanishing surface tension and 
velocity. The quantity in (4.1) is the ^-dependent ki- 
netic coefficient In the notation of Ref. 4, 
6 m =oA <l-cosi»0). From here on, we shall specialize 
to the case m =4, and write *4"=8ftcos 2 0sin 2 0. 

The Ivanisov limit of (4.1), * o (0)~cos 3 0 is a parabola 
which has unit radius of curvature at its tip and which 
trivially satisfies the needle-crystal conditions (2.2) for ail 
v. In analogy to (2.7), let 

*-cqs 3 0=v*, , (4.2) 

and rewrite (4.1) in terms of The result is an equation 
of the form (3.1) with 

PoM 5 tan*+ , (4.3) 



fo(9). 



1 

cos 7 * 



7* * 



r 0 W= 3cos0(5cos*-4) 



8ft 



cos 1 * d$ 

Inserting these functions into (3. 12), we obtain 

XexpnOftstoW*^, 
where A. = -~ + 1 2ft and 



cos" 2 * 



(4.4) 



(4.5) 



(4.6) 
> 

(4.7) 



Notice that, in contrast to the geometrical model (3.16), 0 
in (4.7) is not quite the same as the unperturbed arc length 
£ defined by 



(4.8) 



The integral on the right-hand side of (4.6) can be evaluat- 
ed asymptotically in the limit of small v. Look, for the 
moment, at the case ft =0, so that the only integrals need- 
ed are Fourier transforms of powers of the function 
cosffty). An analysis which is outlined in the Appendix 
yields the estimate 

f* m rf#cos*)V*^C(M>v-°exp( -tf 4 /v^v) , (4.9) 
where 



2ir 
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(4.12) 



7 2 



(4. ID) 
(4.11) 



and 



T(2^/7) ' 

As in (3.16), the dominant exponential in (4.9) is in- 
dependent of fi, but the power of v" 1 in the prefactor in. 
creases with increasing ft. For ft«0, the leading u is ^ 
and * ' 



B 

y29/2B 



■exp - 



I ^ r 



with 



*»-£c<-£)~0.4703, 



(4.13) 



(4.14) 



Both numerical evaluation of the integral in (4.6) and nu- 
merical integration of the fufly nonlinear differential 
equation (4.1) confirm the general form of (4J3). Specifi- 
cally, a graph of 



-^v 3 ^— m [^(0)J^a-2aV^ + . . . 



(4.15) 



versus V^, evaluated for 0.0015 £V£0.01 and extrapolat- 
ed to small Vv, yields a =0. 71 ±0.02=*a, and a value of a 
of about unity, consistent with A. value of B of about 
0.8 fits the integration of (4.1), whereas numerical evalua- 
tion of the approximation (4.6) confirms (4.14) in giving a 
value of B of about 0.5. The latter discrepancy is in the 
same sense and of about the same size as that which oc 
curred for the geometrical model. 

The idea that the dominant behavior of *'(0) at small v 
is determined by the highest power of cos* in (4.6) pro- 
vides some interesting information about the effect of 
crystalline anisotropy. When the anisotropy coefficient ft 
is nonzero, the dominant B dependence in r o (0) comes 
from the second term on the right-hand side of (4.5) and 
enters with the opposite sign. That is, the leading term in 
r 0 {6) is proportional to -ftcos 5 0 instead of +cos 3 0. As 
a result, ir*{0» must be positive as v approaches zero. At 
larger values of v, however, the exponential factor in (4.6) 
is no longer oscillating so rapidly* the slower variations in 
r 0 {6) make the dominant contributions, and <r*(0) be- 
comes negative again. The situation is illustrated in Fig. 2 
where *W) is drawn as a function of v for ft=0, 0. 1, and 
0.2. These graphs extend out to v=0. 1, well beyond the 
asymptotic regime, and have been obtained by integrating 
the fully nonlinear equation (4.1). Equation (4.6) does 
seem to remain a qualitatively good approximation as 
shown by the dashed curves in the figure. 

The important point is that, for ft^O, jr*{0) vanishes at 
a finite value of v, indicating the existence of a steady- 
state needle crystal. For this particular version of the 
model, that is, the minimal boundary*layer model (A— »0) 
with the special form of kinetic anisotropy indicated in 
(4.1), such steady-state solutions occur at arbitrarily small 
but nonzero values of the anisotropy coefficient ft. How- 
ever, the selected value of v=</ 0 u/2>A 5 becomes small as 
ft decreases, and it seems highly unlikely that needle crys- 
tals with very small v or, equivalently, small d 0 can be 
stable. Thus, it seems that the minimum value of ft re- 
quired for dendritic behavior in this model is nonzero and 
is determined by a stability requirement. 
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FIG. 2. r'(0) as a function of v for the boundary-iayer model 
with various values of the anisotropy parameter ft. The solid 
curves have been computed by numerical integration of the non- 
linear differential equation (4.1), and the dashed curves by nu- 
merical evaluation of the approximation (4.6). The heavy dots 
indicate values of v at which the solvability condition is satisfied 
for the corresponding values of 0 4 . 



APPENDIX: SOME MATHEMATICAL DETAILS 

Asymptotic estimates for *'(0). To derive the asymptot- 
ic relation (3,16) for the geometrical model, we must 
evaluate integrals of the form 



<A1) 



in the limit v-*0. This integration is most conveniently 
carried out by transforming to i7=tan0, where 0 is the 
original angular variable that satisfies cos0=sech£. This 
technique is slightly more cumbersome than necessary for 
the geometrical model, but turns out to be specially suited 
for both the boundary-layer model to be considered next 
and the fully nonlocal problem. 1 1 In terms of if, we have 



exp 



where 



(A2) 



(A3) 



Inspection of (A2) and (A3) indicates that we should de- 
form the contour of integration so as to pass near 77=/ 
and then carry out a steepest-descent calculation in the 
neighborhood of that point. (The procedure in this case 
requires that the path of steepest-descent pass through the 
branch cut and back out again,) If |a>| «r 

then 
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0 U+i7?) ,/a J ° (2/fi>') ,/J 
«-f+(2to>»". < A4) 



and 



1/2 



2ir 



exp - 



2Vv 



(A5) 



Note that, as mentioned previously, the dominant power 
of v in the pref actor comes from the largest value of ir 
thus, in (3.16), F * 



#c'(0)«-/ 7/1 <v) 



<A6) 



The same technique works for the integral in (4.9) 
which we shall denote by F^M. We have . ' 

F,(v)= £ m ^_I__ eip |^ (f 1 J , 



(A7) 



with 



*mW- j' d * { Jp?n = fidvtlWf" • CAB 

Again, writing 17 and integrating along the path of 
steepest descent through 77— f, we find 



I 



Xexp 



7Vv 



(A9) 



which can be evaluated without further approximation to 
obtain the results shown in Eqs. (4.9)— (4.12). 

Numerical methods.. Numerical integrations of the 
nonlinear differential equations (2.5) and (4.1) were car- 
ried out with an implicit scheme using 6, *, and 
k=*d*/d§ as functions of arc length £. Initial conditions 
were computed using the asymptotic expansions as close 
to 0=ir/2 as possible, and then the equations were in- 
tegrated back to 0=0. Because, in practice, one can never 
start precisely on the trajectory which enters the fixed 
point, there will always be some components of the oscil- 
lating homogeneous solutions (33) in the function being 
computed. The frequency of the oscillation in (3.3) pro- 
vides an estimate of the minimum step size d£ required to 
resolve these oscillations accurately. In addition, the rate 
at which the magnitudes of these homogeneous solutions 
decay in going toward 0=0 from provides an esti- 

mate of how accurately one must locate the correct trajec- 
tory near tr/2 in order to obtain a desired accuracy of the 
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solution at 0=0. These estimates were confirmed numeri- 
cally (roughly) and were useful for achieving convergence 
of the numerical procedure. 

Numerical evaluations of the various forms of the sol- 
vability formula (3.12), specifically (3.16) and (4.6), were 
performed using (he trapezoidal rule with uniformly 
spaced intervals in (as If performing a numerical 
Fourier transform), . Results were extrapolated to 
dtp/Yv-+Q and tested for convergence in the outer cut- 
off. 

Note added in proof. The connection between singular 
perturbations and solvability conditions has been dis- 
cussed by Barenblatt and Zel'dovich 1 *" in the general 
context of similarity solutions of partial differential equa- 
tions. A recent development along these lines, closely re- 
lated to the work described in the present paper, is the 
discovery of a solvability condition for pattern selection in 
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the theory of viscous fingering. An excellent review of 
the latter topic has been prepared by Bensimon, KadanofT 
Liang, Shraiman, and Tang. 14 I am grateful to L. KanJ 
danoff for informing me about Reft. 12 and 13 and for 
sending me a copy of Ref. 14 prior to publication. 
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Development and evaluation of a multifunctional directly 
compressible diluent consisting of brittle and ductile materials. 
Pharmaceutical Technology , Dec, 2003, by Mukesh C f gQtel, Mpcrnan P, Patei, ShltaLHL 
Bariva. Rlkita K. Dave , Nehaj H, Barjyq 

The authors developed a coprocessed adjuvant consisting of 
mlcrocrystalllne cellulose, dibasic calcium phosphate dihydrate, and 
croscarmellose sodium. The properties of each component were 
evaluated and tablets of nfmesullde were prepared to evaluate the 
functionalities of the coprocessed diluent. It was determined that the 
developed adjuvant exhibited satisfactory tableting characteristics. 

********** 

Direct compression Is widely used In tableting because it requires 
fewer processing steps, Is simpler to validate, and improves drug 
stability when compared with the wet granulation method (1). Direct 
compression also eliminates exposure to heat and moisture during 
processing and is a more economical process. However, the majority 
of active pharmaceutical Ingredients exhibit poor compressibility. 
Therefore, the addition of directly compressible adjuvants Is 
mandatory in such cases. 

A directly compressible filler-binder must exhibit good flowability and 
compatibility. Good flowability Is necessary to ensure rapid and 
uniform die filling, whereas high compatibility is necessary to 
produce tablets having sufficient mechanical strength (2). No single 
diluent is likely to possess alt the Ideal characteristics. For this reason, 
the current trend in Industry is to use multifunctional coprocessed 
adjuvants. 

A survey conducted by Shangraw and Derma rest revealed that the 
most commonly used fillers for tablets and capsules are lactose, 
microcrystalline cellulose (MCC), starch, and dibasic calcium 
phosphate dihydrate (DCP) (3). MCC is well known as a tablet diluent, 
binder, and dfslntegrant. It is a ductile material that shows plastic 
deformation on compaction, which helps produce stronger tablets. 
However, MCC also Is sensitive to hydrophobic lubricants and shows 
poor flowability, and therefore must be modified into large particles or 
agglomerates to Improve its flow characteristics. 
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DCP has a high fragmentation propensity and therefore is insensitive 
to hydrophobic lubricants. It Is also heat stable and has good flow 
properties (4). Other advantages of DCP are its low hygroscopidty 
and low cost. Sangekar et at. evaluated the performance of eight 
directly compressible exclplents. Based on the results of hardness, 
change in disintegration time, moisture uptake, and change In 
volume, DCP was reported as the preferred directly compressible 
adjuvant (5). DCP is best used in direct compression when combined 
with MCC (6). 

Several coprocessed filler-binders are commercially available, 
including Cei lactose ([alphaj-lactose monohydrate and powdered 
cellulose 75:25), Mlcrocelac ([atpha]-lactose monohydrate and 
powdered cellulose 75:25), Ludipress (93% {alpha]-lactose 
monohydrate, 3.5% polyvinylpyrrolidone, and 3.5% crospovidone), 
and Pharmatose DCL 40 (95% anhydrous [beta]-lactose and 5% 
lactitol) (7). The main objective of developing a coprocessed adjuvant 
is to obtain the advantages of each component in a single adjuvant. A 
combination of MCC, DCP, and croscarmellose sodium has not been 
fully explored either by industry or by researchers as a directly 
compressible multifunctional diluent. Therefore, in this investigation a 
coprocessed adjuvant consisting of the two diluents and 
croscarmellose sodium was developed. 

Nearly all excipients for direct compression are manufactured by 
granulation, agglomeration, cocrystaillzation, or spray drying. Spray 
drying has been used commercially to manufacture excipients with 
good Mowability and compressibility, but spray-dried products exhibit 
poor re workability. Crystallization requires critical control of 
processing parameters. The method of wet granulation was adopted 
in this study because of its simplicity. The efficacy of the adjuvant 
was assessed by developing quick-disintegrating tablets of 
nimesulide. 

Materials and methods 

Materials. Nimesulide (100#) was received as a gift from Redson 
Pharmaceuticals Ltd. (Ahmedabad, India). MCC IP (MCC, 100#), 
starch, polyvinylpyrrolidone (PVP K30), hydroxypropyl methylcellutose 
(HPMC K15), and polyethylene glycol (PEG 4000) were received as 
gifts from Zydus Cadila Health Care Ltd.(Ahmedabad). DCP IP (DCP, 
100#) was provided free of charge by Enar Chemle Pvt. Ltd. (Navsarf, 
India). Croscarmellose sodium, sodium starch- gfycolate, and 
crospovidone were received as gift samples from Marutf Chemicals 
(Ahmedabad). Cab-O-Sfi M5 was a gift from Cabot Sanmar Ltd. 
(Madras, India). Magnesium stearate (Vikas Pharma, Ahmedabad) 
and talc (K Chemicals, Ahmedabad) were of IP grade. 

Methods, Evaluation of binders. The premlx of MCC and DCP (70:30, 
50 g) was sieved and then remixed. The blend was granulated using 
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starch paste (Batch [A.sub.l], 10% w/v), PEG 4000 (Batch [A.sub.2], 
26% w/w PEG 4000), PVP (Batch [A.sub.3, 10% w/v PVP), or HPMC 
(Batch [A.sub.4], 0.5% w/v HPMC) as a binder* A melt granulation 
technique was adopted for PEG 4000 and a wet-granulation technique 
was used for the rest of the binders. The liquid binder was mixed with 
the powder blend for 3 mln and then the damp mass was warmed at 
60 [degrees]C for 15-30 mln in a hot-air oven to facilitate 
granulation. The wet coherent mass was then passed through a 44 
mesh. The granules were dried at 60 [degreesJC for 90 mln in a hot- 
air oven. The dried granules were again passed through a 44 screen. 
Fines were removed by shifting the granules on a 120 mesh. The 
granules were evaluated for angle of repose and Carr's Index. Talc 
(2%) and magnesium stearate (1%) were mixed with the granules, 
and tablets with an average weight of 220 mg were prepared on a 
single-punch tablet machine (Cadmach Machinery Ltd., Ahmedabad). 
Tablets were evaluated for crushing strength and disintegration and 
the results are shown in Table I. 

Angle of repose. Angle of repose was measured using the fixed-funnel 
method with a reposograph (model 640, Enar Foundation Research 
Center, Navsari) (8). Per the manufacturer's instruction, the samples 
were graded as excellent, good, fair, or meager if the angle of repose 
was found to be in the range of 30-32[degrees], 32- 35 [degrees], 35- 
37[degrees], or 37-45[degrees], respectively. 

Bulk density and tapped density. The granules were filled in a 100-mL ' 
capacity measuring cylinder up to the 75-mL mark. Bulk density is 
the quotient of weight to the volume of a sample. Tapped density is - 
the quotient of weight of a sample to its volume after tapping the 
measuring cylinder 500 times at a height of ~ 1.5 in. 

Caries index. The percentage compressibility (Carr's index) was 
calculated as 100 times the ratio of the difference between tapped 
density and bulk density to the tapped density (9). 

Crashing strength. The crushing strength of the tablets was 
determined after 24 h of compression (time for stress relaxation) 
using a Monsanto hardness tester (Shital Scientific Industries, 
Mumbal, India). 

Friability. Friability was evaluated from percentage weight loss of 20 
tablets tumbled In a friabilator (model EF2, Electrolab, Mumbai) at 25 
rpm for 4 mln. The tablets then were dedusted, and the loss In weight 
caused by fracture or abrasion was recorded as percentage weight 
loss. 

Disintegration time. The time required for disintegration of six tablets 
placed in the tubes of a disintegration test apparatus (model ED2, 
Electrolab, Mumbai) was measured at 37 (+ or -] 2 [degreesJC using 
900 mL distilled water (10). 
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Evaluation of superdisintegrants. CroscarmeJIose sodium (CCS), 
sodium starch glycolate (SSG) and crospovidone (CRP) were tried as 
superdisintegrants at a level of 4%. Granules were prepared using 
starch paste (10% w/v) and were compressed on a single-punch 
tablet machine. Each tablet was evaluated for disintegration in a petri 
dish (diameter = 5.5 an). A tablet was placed carefully in the center 
of a petri dish containing 20 mL of water. The tablet was allowed to 
wet for 30 s and then the petri dish was subjected to 10 oscillations 
of one Inch each. Particle-size measurements were performed using a 
microscopic method after decanting the clear liquid. The results are 
shown in Figure 1. 

Preparation of physical blend of MCC and colloidal silicon dioxide. A 
physical blend of 97% MCC and 3% colloidal silicon dioxide was 
prepared by blending them in a glass mortar for 30 min. The silica- 
treated MCC (35 g), DCP (15 g), and croscarmellose sodium (2 g) 
were blended subsequently. The granules of 44/120 # were prepared 
using starch paste (10% w/v) as a binder and then were evaluated 
for angle of repose and Carr's Index. Tablets containing silica-treated 
MCC and untreated MCC were separately prepared and compared 
(see Table II, Batch B). 

Lubricant sensitivity test. The effect of magnesium stearate on the 
crushing strength and the disintegration of tablets containing treated 
MCC as well as untreated MCC were studied (batches [C.sub.l]- 
[C.sub.4]). Magnesium stearate (1%) was mixed with the granules 
for 2 min in a glass jar. The granules ready for compression were 
compressed Into tablets and the tablets were evaluated for crushing 
strength and disintegration. Table III shows the results of the 
lubricant sensitivity test. 
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